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Abstract 

We show in detail how the presence of a heat bath of photons effectively gives charged 
particles in the final state of a decay process a temperature-dependent mass, and 
changes the effective strength of the force responsible for the decay. At low tempera- 
ture, gauge invariance causes both these effects to be largely cancelled by absorption of 
photons from the heat bath and by stimulated emission into it, but at high temperature 
the temperature-dependent mass is the dominant feature. 



CERN-TH-6384/92 
February 1992 



1 Introduction 

In thermal field theory, one calculates thermal averages of observables Q: 

(Q) = Z- 1 Y,( i \^ H Q\i) (!) 

where the sum extends over a complete set of physical states. As usual, H, Z and (3 are respectively the 
Hamiltonian, the grand partition function and the inverse temperature of the system. The perturbation 
theory that can be developed to calculate such quantitities (Q) is very similar to the familiar Feynman 
perturbation theory W . Thus it seems natural to suppose that concepts from Feynman perturbation 
theory apply also to calculations in the presence of a heat bath. A particular example is the idea of 
mass renormalisation and the supposition that thermal effects somehow give a particle an additional 
"thermal" mass, though there has been disagreement about how exactly this additional mass should 
be defined PL In this paper, we derive the mass correction, and consider the other effects that the heat 
bath has on reaction rates. 

Since our aim is to elucidate basic principles and methods, we study a very simple problem. We 
consider a heat bath of photons in equilibrium at temperature T and a neutral scalar particle $ that 
decayst into a charged scalar particle and its antiparticle. We take the decay interaction to be through 
the weak pointlike coupling \$<ft <j> and work to lowest order in A. For simplicity, initially we do not 
work with proper scalar electrodynamics, but pretend that the photons are scalar particles with a 
mass jj, and the simple interaction 

(eA + e 2 A 2 )<pU 
Later, we go over to the proper gauge theory. 

We suppose that the heat bath has not been in existence for a time long enough for the particles <!> 
and (f> to be anywhere near thermal equilibrium; rather, the heat bath will be supposed to consist only 
of the photons. This is again because we want to consider the simplest possible problem in order to 
concentrate on basic principles. As we describe briefly at the end of section 3, it is straightforward to 
extend the analysis so as to include the particles <1? and <f> in thermal equilibrium in the heat bath. 

In section 2 we show in detail how the heat bath effectively gives the charged particles an additional 
squared mass of order e 2 T 2 , and how it also changes the effective value of the decay coupling constant 
A. Further, there are effects arising from absorption of photons from the heat bath, and stimulated 
emission of photons into it. 

When we go over to the proper gauge theory in section 3, additional features enter. This is that gauge 
invariance ensures that there are various cancellations. Part of the cancellation is of the infrared 
divergences that are now present in the change in the effective value of A and in the absorption 
and stimulated emission! 5 !. But we find that it is not only the divergences that cancel: in the low- 
temperature expansion up to order e 2 T 2 everything cancels, and the net effect of the heat bath on the 
decay rate is zero up to this order. However, at high temperature the mass-shift effect is of order e 2 T 2 , 
while the other effects together are only of order e 2 T, so that the mass shift is the most important 
feature. 

2 Change of effective mass and coupling 

We calculate to order e 2 the difference between the decay rate in the heat bath and that at zero 
temperature (which is the same as that with no heat bath, as calculated in ordinary Feynman theory). 
It is convenient to work with real-time thermal perturbation theory for the electromagnetic interac- 
tion. Normally, in real-time thermal perturbation theory there are complications associated with field 
doubling! 3 !, but these are not encountered when the heat bath contains only photons. 

The derivation of the real-time perturbation theory requires an adiabatic switching off of the interaction 
with the heat bath at asymptotic times^. We impose the switching off in the rest frame of the heat 

t Our results apply equally well when the decaying particle has spin. As is customary, we shall ignore the 
problem that strictly-speaking there is no proper field theory of unstable particles. We could equally 



bath (though the choice of frame for this probably is not essential). We do not assume that the 
decaying particle is necessarily at rest in this frame. We shall find it necessary to retain explicitly 
the switching factor in the first part of our calculation. We shall take it to be S v (t). We also need 
to introduce a switching factor S^(t) for the interaction that causes the decay of the $ particle. Our 
results will not be sensitive to the precise shape of these functions, provided they are very close to 
unity for all finite t and go to zero when \t\ — > oo. It will be helpful to have in mind the explicit forms 
Sri(t) = e ^1*1 and S^(t) = e~^, with Fourier transforms 



Eventually, we let 77, £ — > 0, and then both functions in (2) become 2it5(uj). But before we take that 
limit, at each vertex of our Feynman graphs there will be an additional incoming energy to, with a 
corresponding multiplicative factor (2) and an integration over oj. The heat bath is coupled to the 
decay process for a very long time when 77 is very small, but as it is nevertheless switched off in 
the remote past and future the mass of the (f> particle in the asymptotic in and out states is just its 
zero-temperature value m. We show that nevertheless the self-energy insertions generated by the heat 
bath effectively change this value in the formula for the decay rate, and derive the other effects caused 
by the heat bath. 

A change of mass is, of course, usually associated with renormalisation. We recall that the familiar 
zero-temperature renormalisation is partly a mere rearrangement of the unrenormalised perturbation 
expansion, but also it involves making a change in it: one introduces a new choice of the asymptotic 
condition to be satisfied by the fields. The mere rearrangement concerns how one deals with the self- 
energy insertions on internal-line propagators in an amplitude, and has two aspects. A part 5m 2 of the 
self-energy H(p 2 ) is subtracted from II and added to the squared mass: m 2 — ► m R = m 2 + 5m 2 ; this 
subtraction is chosen such that the new II vanishes at p 2 = m 2 R , so that the propagator pole is at this 
point. Also, Z = (1— lT(m|j.)) _1 is factored off from the propagator and the square root of this factor is 
instead made to multiply the coupling constant at each end of the propagator. If all the vertices of the 
Feynman graph were internal vertices, this mass and coupling-constant renormalisation would be just 
a rearrangement that does not change the Feynman integral. But in order to justify the same change of 
coupling constant at external vertices, where there are a smaller number of propagator lines attached 
than at internal vertices, the asymptotic condition must be changed so that the normalisation of the in 
and out free fields is changed by a factor y/Z. Also, these fields are required to have mass m R instead 
of the original mass. This change of asymptotic condition, being more than just a rearrangement, must 
be given physical justification; in electrodynamics one usually does this by noting that the lowest-order 
electron-electron scattering amplitude at small momentum transfer, renormalised in this way, agrees 
with the Rutherford formula, say, when one identifies the renormalised mass and charge with those 
listed in the data tables. Happily, the renormalisation absorbs the infinities of the theory into the 
renormalised quantities, but one would need to go through such a procedure even if there were no 
infinities. There are many other renormalisation procedures that absorb the infinities equally well, 
but the one we have described is the one that most directly makes contact with the data-table values 
of the physical quantities. 

There is obviously no reason why one should not make an exactly similar rearrangement in thermal 
perturbation theory, where now the self-energy has an additional temperature-dependent part, since it 
does not change anything. That is, one is free to rearrange the calculation such that a new temperature- 
dependent mass and a new temperature-dependent coupling are used on internal lines and vertices. 
This mass will correspond to the position of the pole in external legs of the Green's functions of the 
finite-temperature theory, and the coupling will correspond to a pole residue. But some justification 
is needed if one wants to assume that these new quantities enter in a simple way in the calculation of 
finite-temperature reaction rates. We shall show in this paper that in fact one does not need to make 
any new assumptions to justify this. 





(2) 



We set up the perturbation theory in the interaction picture that coincides with the Heisenberg picture 

AHV. 



with the particle initially, it is nevertheless present and so we should describe the initial system by 
the density matrix 

p = Z - 1 ^e-^at(P) a*(P) (3) 

i 

where \i) are the possible states of the heat bath of noninteracting photons and Hq and Zq are its 
Hamiltonian and grand partition function. The decay probability is 

£</|TpTt|/) = Z-Hr e-^o a$ (P)rrt4(P) (4) 

/ 

where we have used the completeness relation for the states |/) and the trace is in the subspace of 
heat-bath states. We insert the standard perturbation-theory operator expression for the transition 
matrix T, using the in fields, and evaluate the terms of (4) that are of order e 2 A 2 . This results in the 
types of diagrams of figure 1. In these diagrams the charged-particle lines correspond to the usual 
zero-temperature propagators and 5-functions, with the zero-temperature data-table masses, because 
the heat bath is supposed not to contain any charged particles. The photon lines correspond to 

A(k;(3) = J d A x e- ik - x Z~hv e~ pH S A(x)A(0) (5a) 

for the real-photon lines and 

D(k;0) = J d 4 x e- ik - x Z-Hr e~^ H o T (A(x)A(O)) (56) 

for the virtual ones (with T the real-time time-ordering operator). (Recall that, initially, we are taking 
the "photons" to be scalar particles with mass fj,.) 

Because we want to calculate the difference between the finite-temperature and zero-temperature 
decay rates, we subtract from (5) the zero-temperature expressions, namely 

A(fc; oo) = 9(k )2ir5(k 2 - /i 2 ) (6a) 

and 

= < 66 > 

After this subtraction, we are left with 

A?{k) = D?(k) = 2ir5(k 2 - fi 2 )fP(k ) (7) 

with 

fP(h„\ - 

oP\k 



f"(ko) = r ( 8 ) 



Although these two functions are equal, they represent different physical processes: A@ corresponds 
to permanent absorption or emission of a photon from or into the heat bath, while D 13 represents 
absorption together with emission (or vice versa). Notice that, unlike at zero temperature, the real- 
photon function does not contain a factor 9(k ); this is because it describes both emission of real 
photons into the heat bath and absorption from it, depending on the sign of k . 

The additions to the two self-energies appearing in figures la and lb arising from the heat bath are 

n f ( Pl ) = e 2 f -^-DP(k)- -i 5 (9) 

byF ' J (2vr) 4 v J ( Pl + k) 2 -m 2 + ie y) 



At each vertex of the graphs of figure 1 there is conservation of 3-momentum. However, we must 
allow for a spurion energy u, associated with the switching factors (2). Although eventually 77, £ — > 
so that the spurions then disappear, we need them initially to make sense of the first two graphs. 
Without the spurion energies, the internal (zero-temperature) scalar-particle propagator would be 
(p 2 — to 2 ) -1 which then has to be evaluated at its pole because there is also an on-shell line with the 
same momentum. But before 77, £ — > the on-shell line's momentum is not the same as that of the 
internal-line propagator and so the pole contribution is regularised. It will be a useful check on our 
final answer that it is independent of the relative rates at which 77, £ — ► 0, or indeed of the precise 
shapes of the switching-off functions. 

Each graph of figure 1 represents an amplitude on the left-hand side of the cut lines and the complex 
conjugate of an amplitude on the right. Denote the spurion energy leaving the left-hand decay vertex 
by Q and that entering the right-hand one by $7'. Similarly, introduce energies u for the spurion 
energies that leave the heat-bath- interaction vertices to the left of the cut lines, and energies to' for 
those that enter those to the right. The squared matrix element contains, among other things, 

/ § § n (£) n (£) ** (° + em - (« + em - 



s £ (a)S!(n')n s i( u >n s >') 

where the switching-off functions S are defined in (2) and AP is the difference between the initial 
energy and the sum of the energies of the cut lines. 

When we calculate a probability from a squared T-matrix element, we normally have a factor 
(27r) 4 <5( 4 ) (AP) in both the amplitude and the complex-conjugate amplitude, where AP is the difference 
between the total 4-momentum on the cut lines and that on the external lines. In order to obtain a 
cross-section we normally factor off (2tt) 4 5 ( -^ (0), which is interpreted as the volume times the total 
interaction time. This last has been affected by our introduction of the switching factors: while we 
still have (27r) 3 <5( 3 ) (AP) in both the amplitude and the complex-conjugate amplitude, the square of 
2tt5(AE) is no longer there. In order to see what it has become, consider a graph that does not "need" 
the switching factors, that is one for which we may set the oj and oj' to zero in the propagator factors 
before performing the integrations in (10). The integral (10) is then equal to 

J dtdt'e^^S^S^Hs^Hs^t') (11) 

In the limit where the switching factors and S$ become equal to unity this is just 2tt5(AE)2tt5(0). 
Before we take this limit, (11) is small except very near AE = 0, that is it is almost equal to 2ir5(AE) 
times a large factor N obtained by integrating (11) with respect to AP: 

N n = J dt[S^t)] 2 [S v (t)] n (12a) 

where n is the total number of heat-bath vertices. It is this integral that we must identify as the 
total interaction time and factor it off. In the case of the simple explicit forms (2) for the switching 
functions it is 

2£ + 7177 

where we have omitted terms that remain finite or vanish when £, 77 — > 0. 

Consider now the graphs of figures la and lb. When we insert the spurion energies, the total energy 
Ei + P 2 on the cut lines is not equal to the initial energy P, because there is no delta function 
8{E\ + P2 — P). But it is convenient to insert the delta function 8{E\ + P2 + AP — P) and remove it 
again by integrating over AP. Since the switching-off functions constrain the spurion energies to be 
small, we may expand this delta function in a Taylor series: 



5(E 1 + P 2 + AP - P) = — f {AE) V 8^ (E 1 + P 2 - P) 



Then the expressions for figures la and lb give contributions to the decay rate that are of the form 
S /^2rf(+)(p^m 2 )^2^+Hp^ m 2 )(2»)'S«'(pi+P2-P)E Ml '' 1 ( £ l+^-- E )'"(Pl) 

For figure la 
and for figure lb 

^W ^^ji^ ^^ri (146) 

We must add on the complex conjugate of each of these, corresponding to the self-energy insertion 
being made instead in the complex conjugate amplitude. With the explicit forms (2) the integrations 
in (14) are easy to perform by taking residues at the poles of the integrand. Then, when 5,)j^0we 
find that for both figure la and figure lb we may write the result as 

W» fn ^ - n/ V) , d n/ V) 



2E{ dEf 
IIP fa) 



2E l 

iff fa) =0 for v > 1 (15) 

(though note that for figure la IP is independent of p\ and so the derivative term vanishes). In the 
Appendix, we show that the results (15) are general in that they do not require the explicit form (2) 
for the switching functions. 

(13) and (15) correspond to self-energy insertions in the line p\\ there are similar terms for insertions 
in the line p 2 . 

The derivative term in Iq is naturally interpreted as the finite-temperature renormalisation of the 
squared coupling A 2 by the self-energy insertion in figure lb. At zero temperature the self-energy is 
a function of pi and the coupling renormalisation involves its derivative, which is naturally regarded 
as being with respect to pi, but one could equally well use E\ because the result is the same. At 
finite temperature the heat bath breaks the Lorentz invariance and the variable of differentiation is 
necessarily E\. The graph of figure lc also may be thought of as just further contributing to the 
change of the coupling A. 

The other terms in (15) shift the mass m. To see this, note that the zero-temperature decay rate T 
is 

r ° = S / |^|^ (27r)V4)(pi +P2 - p ^ 6i+) & - ™ 2 )^ {+ HpI - - 2 ) 

A 2 f d 3 Pl 1 



Because 



2E J (2vr) 3 AE x E 2 
d 1 1/1 1 



IUf utE^ S{El+E >- E) <16) 



dm 2 E X E 2 ~ 2E 1 E 2 \E 2 E\ / 
JL S{El + E 2 -E)= 5'(E 1 + E 2 -E) f -i- + -M 



when we insert into (13) the nonderivative term of l(j and /f from (15), together with the corresponding 
contributions from the self-energy insertion instead being made in p 2 , we obtain a correction to the 
decay rate equal to 



Ar (; 



— j 

2E J 



d 4 pi d 4 p 2 
(2vr) 4 (2vr) 4 



(27r)V 4 )(p 1+P2 -P) 



where 







*m 2 (pi)^27rf<+>(p?- 



m 



2tt5^ (pi - m 2 ) + (px < — ► p 2 ) \ (17a) 



5m 2 ( P )= nf(p) + nf(p) 



(176) 



Note that, while the contribution to (5m 2 from nf (p) depends only on (3, that from IT^(p) depends 
also on |p|. The effect of the heat bath on the decay rate depends on how fast the decaying particle 
is moving. 

We have shown, then, that the first three graphs of figure 1 may be thought of as changing the mass of 
the decay-product particles and the decay coupling. This has not needed any new assumptions about 
how to renormalise at finite temperature, but is a consequence of standard zero-temperature field 
theory combined with statistical mechanics. In addition, there are the graphs of figures Id and le, 
which each represent the absorption of a real photon from the heat bath and the stimulated emission 
of a real photon into it. 



3 Scalar electrodynamics 

Consider now the proper scalar electrodynamics, with zero-mass photons. It is convenient to use the 
gauge Aq = 0. Then the thermal parts (7) of the real and virtual photon propagators become 



Af#) = flg(fc) = (<% - ^) 2n5(k 2 )f(k ) 



They are purely transverse: in the real-time formalism the longitudinal-photon propagator does not 
have a thermal part^. The self-energies (9) become 



nf(Pi) = e 2 / I 0^(fc) = ie 2 T 2 



where we have used 



n? to) - / (gi 4 Pl .Df # K (pi+t) , 1 _ m , + fe (IB) 

x dx 



f 

Jo 



e x - 1 6 



2 



and written 1/(3 = T. Again, nf is independent of p\. 

The Aq = gauge is particularly convenient because with it the spurion energies do not appear in 
the numerator factors of the Feynman integrands, only in the denominators. So the analysis of figures 
la and lb is similar to before. However, there are some differences. First, when we set p\ = m 2 as 
required by the (5-function in (13), we obtain 

Df(Pi) = e ' / ^uD^{k)p l3 {p -L^ - m5{2 Pl .k)} (19) 

The integrand for Re 11^ is odd in k and so the integral vanishes. 11^ also has an imaginary part, 
which arises because the Bose distribution contains an infinite number of zero-energy photons; it is 
proportional to T. (In the previous section, where we pretended that the photon had a nonzero mass, 
this imaginary part was zero.) However, the imaginary part cancels when we add to figure lb its 

t-\1 cs-v /~t\i~i tit rrofo Time -tTT-Vi<-\1/zi /~i"f -j-Vin nrlrlif lnnnl wiocc oVinf-l- ( 1 TVv\ nnnonrl Vvt t It q a f V\o-|-V* ^nvnflo 



from the seagull self-energy insertion, figure la. This is peculiar to scalar electrodynamics in the 
A = gauge. Gauge invariance ensures that we obtain the same total thermal mass shift for scalar 
electrodynamics in other gauges, but it generally does not come only from figure la: in other gauges 
(p) does not vanish at p 2 = m 2 . Also the spurion energies may appear in the numerators of the 
various integrals. The same applies to spinor electrodynamics, even in A = gauge. 

The other point is that the derivative term in I^{pi) in (15) is now infrared divergent when we set 
p\ = m 2 as required by the 5-function in (13). This infrared divergence is cancelled by the graph 
of figure Id. In order to handle the infrared divergences, we add together the graphs before we take 
the limits ^, 77 — > 0. There are also infrared divergences in figures lc and le; again they cancel each 
other. In each of the two cancellations a finite piece is left over. This finite piece is a combination of 
the change in the effective value of the decay coupling constant A and the absorption and stimulated 
emission of real photons. 

Because the mass shift 5m 2 of (17b) arises just from Ilf , which according to (18) is independent of 
momentum, the correction (17a) to the zero-temperature decay rate T is just 5m 2 -§^z- The zero- 
temperature decay rate is 

A 2 



8irME 

and so the mass shift changes this by an amount 



y/\M 2 -m 2 (20) 



2 2 dT° \ 2 e 2 T 2 1 

a ~ ^ dm? ~ 96itME ^^M 2 - m 2 ( ' 



The change in the effective value of A changes the rate by 
A 2 e 2 



AIV = -^7^5 J d 4 kf(k )6(k 2 )d 4 Pl 6M(p 2 - m 2 )d 4 p^ + \p 2 - m 2 )5^(p 1+ p 2 - P) 



- p 2 - ( Pl .k) 2 /k 2 | p 2 - ( P2 .k) 2 /k 2 2 Pl . P2 - Pl .k P2 .k/k 2 

(Pi-k) 2 (pi-k) 2 " (pi-k)(p 2 .k) 

where f@ is the Bose distribution function (8). In the square bracket, the first two terms correspond 
to the self energy insertion of figure lb and similar graphs, while the third corresponds to figure lc. 
The way in which the contributions from these figures distribute themselves is gauge-dependent and 
the separate terms break Lorentz symmetry not only because of the presence of the Bose distribution 
associated with the heat bath, but also because we have used the Aq = gauge. However, if we use 
the 5-functions in (22) we find that we can reduce the square bracket to the Lorentz scalar 



(22) 



_ ( _ P2 \ 

\ Pl .k p 2 .k) 



so that in the sum of the three terms the Lorentz invariance is broken only by the heat bath. It is 
straightforward to verify that the same answer is obtained for Ar& c in Feynman gauge. 

The contributions Ar^ e from figures Id and le, and similar figures, sum to an integral exactly the 
same as (22), except that it has the opposite sign and 

5 (4) (pi + p 2 - P) —> 5 (4) (pi + p 2 + k - P) (24) 



It cancels the infrared divergence in (22). In order to handle this new 5^ -function, it is convenient 
to make changes in the integration variables P \ and p 2 so as to bring it back to its original form. We 
choose 

p 2 -k 

Pl ^ Pl ~Pk k 

Vr> > V-?. zr^k 25 



These changes of variable have unit Jacobian and leave the square bracket (23) unchanged, but they 
result in 



S' + '{pi - m 2 )«<+>(j, 2 - m 2 ) —> (pj - m 2 - 2 2l±?2±\ S (+i _ m ' _ 



(26) 



That is, the difference between and T^e, apart from one of overall sign, is just in their <5-functions 
(26). If we want a low-temperature expansion in the temperature T, we may expand this difference 
as a Taylor series. The first surviving term gives 

A 2 e 2 If d 

Ar 6c + AY de = 



E (2vr) 5 \dm 2 



J d 4 kfP{k )5{k 2 )d* Pl 5^\p\ - m 2 )d 4 p 2 5 {+ \p 2 2 - m 2 )5 (4) ( Pl + p 2 - P) 



M — — — tl m 

(P.k) 2 



X 2 e 2 T 2 1 



96itME yitf 



(27) 



In this approximation, the decay rate in the rest frame of the initial particle does not depend on how 
fast it is moving relative to the heat bath. (That is not true of the contributions from the higher-order 
terms in the expansion of the product of 5-functions in (26), which yield higher powers of T.) 

The result (27) for ATb c + AFde exactly cancels Ar a in (21). That is, the mass-shift effect on the 
decay rate is cancelled, up to terms of order T 2 , by the effects of the change in the effective decay 
coupling and of absorption and stimulated emission, whatever the speed of the initial particle. 



This result is applicable for low temperatures, T <C yj \M 2 — m 2 . If, instead, the temperature is high 
(though still e 2 T < \M 2 — m 2 in order that the perturbation expansion may make sense), we find 
that ATbc + Ar^e is much smaller than Ar a . To see this, note that from (22) and (26) we have 

/k rlk 
^k^jJ(k,P) (28) 

where J is even in k, finite at k = 0, and vanishes as 1/fc 2 for large k. When (3 is small we may 
approximate the denominator in (28) by f3k, so that (28) becomes T times a convergent integral. On 
the other hand*, the result (21) for Ar a is exact for all T. That is, at high temperature the mass shift 
(which in A = gauge arises just from figure la) is the dominant effect, while at low temperature it 
is cancelled by the other effects. 

So far we have supposed that the heat bath contains only photons. If the charged particles also are in 
thermal equilibrium, there is already an effect independent of the interaction with the photons: when, 
as we have been supposing, the charged particles are bosons, there is stimulation of their emission 
into the heat bath. If the initial particle P is at rest relative to the heat bath, this effect multiplies 
Tq by [1 + /^(^M)] 2 . (For fermions, there is a change of sign, so that the heat bath tends to suppress 
the decay.) The corrections in order e 2 are partly zero-temperature renormalisation to the mass and 
decay coupling in this stimulated emission rate T . The analysis in the Appendix of the effect of a 
self-energy insertion applies also when it involves a thermal charged-particle propagator (though now 
one must take account of the field doubling of real-time finite-temperature field theory! 3 !), and so there 
will be additional temperature-dependent corrections to the charged-particle mass and to the decay 
coupling. At low temperature all these effects are exponentially small, so that the 0(e 2 T 4 ) correction 
to T from the photons in the heat bath is the dominant one. 

4 Conclusions 

We have shown that, in order e 2 , the presence of a heat bath induces a mass shift, proportional to 
the square of the temperature, for charged particles that are involved in a reaction in the heat bath. 

* If we had similarly expanded the denominator in the integral that gives Ar a , we should have obtained 

' / ' ^~ i m tie o rln rovrmnf nrrv o 1 



However, this mass shift cannot be dissociated from other effects, including vertex corrections and 
stimulated emission and absorption of real photons. Gauge invariance results in cancellation of the 
infrared divergences that are present in these separate effects. But at low temperature it also leads to 
a cancellation of the mass-shift effect itself, so that the correction to the reaction rate behaves only as 
the fourth power of the temperature. On the other hand, at high temperature the mass shift is the 
dominant effect. 

Our calculations are based on standard Feynman techniques of zero-temperature field theory together 
with statistical mechanics. They are necessarily applied to reaction rates rather than to amplitudes, 
and we find that field theory at finite temperature requires no new assumptions about renormalisation 
or asymptotic conditions. 



We are grateful to Antti Niemi for reading an early version of the manuscript and to Tai Wu for 
reading a later one. 



Appendix 

We derive the result (15) for the integrals (14) for general switching functions S v (t) and S^(t), assuming 
only that they are real and almost equal to unity for all finite t, but go to zero for infinite \t\. 

In (14a), which corresponds to figure la, replace each of the switching functions S by its Fourier- 
integral representation and perform the AE, T' and uj integrations: 

tf(Pi) = ^i" 1 / dTdtS^T) [{-i-^J S € Cn] S v (t) [0(t-T)e-^-^+e(T-t)] Ufa) 

(Al) 

We must add on the complex-conjugate integral. Then, integrating by parts with respect to t and 
using the definition (12a) of Ni we have 

^(Pi) = -2^nf(p!) (A2) 

together with a term proportional to 



/ 



dTdt[S i (T)] 2 ^^-cos(2E 1 (t-T)) (A3) 

This integral is negligibly small: we may take the limit £ — > under the integral because the cosine 
ensures that the T integration converges without this switching factor, and then vanishes. The t 
integration vanishes also. The result (A2) agrees with (15), because for figure la H(pi) is independent 
of pi. By similar manipulations one finds that if (fi) is equal to the result given in (15), together 
again with a negligible term proportional to (A3), and that Ig = for v > 0. 

The corresponding analysis for the integral (14b), corresponding to figure lb, is only a little more 
complicated. Now does depend on p\ . Since for small n, £ only small spurion energies are important, 
we may expand U^(Ei + w,pi) in a Taylor series in powers of u. With similar manipulations, we 
find that only the first two terms in this expansion survive, and we obtain the results in (15). Notice, 
however, that when the photon mass is taken to be zero, as it is in section 3, Jl^(Ei + u, pi) is not 
analytic at u = and so such a Taylor expansion is then suspect. 
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